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Summary, 

On January 3, 1918, the authors presented a paper on *' The Critical 
Loading of Struts and Structures." Since then the method of analysis has 
been extended to include the more general problem of the vibration of a 
continuous beam simply supported at any number of non-collinear points, the 
external loading being a constant longitudinal end thrust, varying from bay 
to bay and a periodic lateral loading. In addition, the supports are assumed 
in a state of periodic vibration. The flexural rigidity constant along each 
bay varies, however, from bay to bay. These conditions will then corre- 
spond closely with those originated in a wing spar of an aeroplane when in 
actual flight and influenced by engine throbbing. A very generalised form 
of the equation of three moments is derived, and the conditions for resonance 
and crippling expressed in a convenient determinantal form. The general 
equation where the restrictions are removed as regards constancy of end 
thrust and of flexural rigidity is then treated and illustrated in the 
particular case of the crippling of a strut of variable flexural rigidity. 



Attention has at intervals been drawn to the effect of vibrations set up in 
aeroplane spars by external forces during flight, such as engine throbbing 
and wind forces. For very simple cases a few investigations exist to find 
under what conditions the results of this become serious on account of 
resonance, But no proper treatment of the actual problem and of the stresses 
brought into being appears to exist, especially in the case where the spar is 
under end thrust. In a previous communication* the static case of a beam 
supported at any number of non-collinear points and under the influence of 
lateral loading and end thrust varying from bay to bay was investigated, and 
the conditions for instability derived in a determinantal form. It will be 
seen that the previous discussion falls out as a mere particular case of the 
more general question developed in the present paper. 

It is clear at the outset that the natural period of vibration and the con- 
ditions of resonance will be considerably affected by the presence of the end 

^ ' Eoy. Soc. Proc.,' A, vol, 94, p. 405 (1918). 
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thrust, but it is not possible from simple considerations to anticipate to what 
extent. In any case, the problem is complicated in aerodynamic practice by 
the fact that the presence of bracing wires attached at intermediate points to 
the spar causes a variation in this thrust along the bay in the vertical plane 
and even in the longitudinal force plane. Under static conditions, it was 
previously understood that a structure of this type would fail through 
instability when the length of any one bay became equal to that of Euler's 
critical length, but it has already been proved in the authors' previous paper 
that this is not the case, and the real conditions under which failure arises 
were there stated. It appeared that the bays could not be regarded as 
isolated struts under end thrust, but that on account of interference effects 
between bays considerable modification arose. It is to be anticipated that in 
the analogous problem here to be investigated, corresponding interference 
due to vibration would result, and that likewise the spar could not be treated 
as a series of isolated struts. under end thrust and external forced oscillation. 
It is proposed for the case of a spar supported simply at any number of 
points to investigate the effect of end thrust differing from bay to bay of the 
spar, upon the natural period of vibration, upon the conditions of resonance, 
and upon the distribution of bending moments and stresses in the structure. 
The ground will be cleared initially, and the manner of treatment exemplified 
by dealing in detail with the case of a beam under end thrust, simply 
supported at both ends, under a vertical periodic lateral loading, and a 
periodic bending moment at one end. 
Let the lateral loading be given by 

w' z=: X COS qt, 

where X is a function of x only, and let the end thrust be F. If Mo and Mi 
be the bending moments at the two ends, and if P be any point in the beam 




of co-ordinates (x, y) at any time t, then the moment on the portion to the 
left of P is 

— M — Mo+Zol/— ^)~-l>— '^|r(a],2;), 

where /o is the shearing force at 0, and i|r(:??, t) is the moment of the lateral 
loading on that portion about the point P. 
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Since the portion of the beam is in motion under the influence of these 
moments, this must be equated to the rate of change of angular momentum 
of that portion, i.e., 

(j ct^ 

where w/g is the mass per unit length of the beam assumed constant along 
the spars and neglecting the rotation of the element. Equating these two 
expressions and differentiating twice with respect to x, we find 

3^M -u.3^?/ d^ __wd^y 

dx^ dx^ dx^ g dt^ ' 

Eemembering that EI ^^-^ = M and that ,~\r is equal to the lateral loading, 
the equation reduces to 

where \^ = F/EI and jjp = w/g'EI. 

The problem will be solved by determining the solution of equation (1) 
which satisfies the conditions— 

(a) ^ — : at ^^ = and x = I for all time, 

(h) 'EI'— = Mo and Mi at ^ = and x ™ / respectively, where Mo and 

Ml may be any given functions of t, 

(c) 2/ = a known function of x for ^ = 0, 
{d) dy/dt = a known function of x for t =: 0. 

Since the equation has constant coefficients, solutions are of coui'se additive. 
Let the expression for y be supposed expressed in the form 

y = l^An cos qj + ^Bn sin qj, (2) 

where An and B„ are functions of x, at present unknown, qn unknown con- 
stants. Inserting (2) in equation (1) it follows that 

s=: X cos qt. (3) 
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It will be sufficient to equate the coefficients of cos gj and sin qj on each 
side of the equation, and we find accordingly 



^"^^-fX^^^y^VAn = 0,"^ 



except where q„ = q, when 



cfo;* da? 



y (4) 



To find a solution of these equations, the end conditions as restrictions on 
y and its derivatives must now be interpreted in terms of the restrictions 
on A and B, etc. Condition (a) which states that y = at (» = and x = I 
for all time requires generally that A« = and B^ = at these points. 

To illustrate the problem it will be sufficient to suppose 

Ml = 'Ela cos pt. 

This corresponds with an applied periodic bending moment at one end, but 
since the solutions of the equations are additive, the presence of one single 
such term will not detract from the generality. Hence 






= at ^ = Z f or all time, 



EI ^ = Ela cos pt Sit x=: for all time. 

ax^ 

These interpreted for the A's and B's require that 

ct A^fij ctOii ^^ u at X — If 

d^Bn/dx^ =: Q at X =: I and a; = 0, 
and d^An/dx^ = at a? = 0, 

except when qn = p, when 

ct J\.pJ CtX '^^ ct. 

When qni^p the solution given by (4) is 

A„ = Pji sinh kn^ + Qn cosh hnX + E„ sin JcnX + S„ cos knX, 
where ±'kn\ ±kni are roots of the equation 

JA^XVc^--fJlPq/=:0, (6) 

so that K' + kn'^ = (XH4/iV)^- 0) 
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Inserting the boundary conditions 
= Qm + S„, 
= (Pn sinh knl + Q„ cosh knl) + (K„ sin kj 4- S„ cos kj), 

= kn^(Pn^inhknl + QnG0S}lknl)-'kn^(BnBinkJ+QnC08kJ), 

from which it follows that P„, Q„, E,„ and S„ are zero, unless 

sin kJ = 0, 

le,, unless kJ = tt, 2'7r, Stt, ... etc., (8) 

where P«, Q«, S^ are then zero and E« is yet to be determined. This last con- 
dition may be written otherwise, viz., 

P 2 

01 g, ____.^__^„^. (9) 

n is any integer, and as long as 

FP ^ 2 2 

EI 

there is a real period. 

For n^l when this becomes an equality, the normal first Euler's critical 
occurs, and above this value of the end thrust the period 27r/qn becomes 
imaginary, so that the value of y increases exponentially, and the motion of 
the beam is unstable. The same argument applies for the higher criticals. 
It should be particularly remarked that the fact that the beam is in a state 
of vibration has not affected the value of the crippling load, but that the 
longitudinal thrust is such as to increase or decrease the period of natural 
vibration according as the force is compressive or tensive. 

In the general expression for y there appear, therefore, only terms in 
cos qj and sin qj for which qn has the value derived in equation (9), 
together with terms in cos pt and sin qt to be obtained as follows : 

So far 

CO . 

2/ == 2 sin knX (Rn cos qj -f E»/ sin qj) + A cos qt + A' gob pt, (10) 



where qn is given by (9). 

The conditions that A and A' are to satisfy are 

(1) At 5^ = and x = I, A and A' are zero, 

d^A .d'A' 

"■"' ■■■' " ^ ' ailtl — z — j: 

dxr ax' 



J 



(2) At X = /, -y-^ and -j-j are zero. 



(3) At a:' = 0, f^ = and ^' = a, 
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where A and A' satisfy the equations 

Consider the expression for A', viz., 

A' = P sinh k'x + Q cosh ¥x + E sin Jcx + S cos fa, 
where ±k' and +M are the roots of the equation 

^.4^.X%2_^2p2^ 0. (13) 

Inserting the end conditions, the following hold : — 
= Q + S, 

= (P sinh l/l + Q cosh ¥1) + (R sin Id + S cos Jd), 
= ¥^ (P sinh ¥l + Q cosh ¥1) - F (E sin M + S cos M). 

from which P, Q, Eand Sare easily determined in general as finite quantities. 
This fixes the term A' cos p^ in (10). Special circumstances arise, however, 
when 

sin M = 0, 
ie., when M = tt, 27r, etc. 

This implies that 2'ir/p, the period of the externally applied bending moment 
at the end synchronises with the natural period of vibration for the whole 
beam under the end thrust. Under these circumstances P, Q, E and S 
become excessively large, and y also as a consequence of the synchronism. 

To determine A we write 

A = Pi sinh kix -{- Qi cQsh kix -f Ei sin kix -f- Si cos kix + 



O 9t > 



where +/<;i' and ±hi is obtained from 

The end conditions are of course A = and d^A/dx^ = for x — and x = I. 
Hence 

= Qi + Si + Z(o), where Z (x) = -—— A 

= Pi sinh h'l + Qi cosh h'l + Ri sin hi + Si cos hi + Z (I), 

= Qih"-Sxh^ + Z"io), 

= h'^ (Pi sinh fci7 + Qi cosh h'l) - h" ( Ri sin hl+ Si cos hi) + Z" (I). 
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From these four equations the values of the four constants Pi, Qi, Ei, Si 
are easily derived as definite finite quantities, except in the special case 

when 

]cil = TT, 27r, etc., 

i.e., when the lateral impressed force synchronises with a natural vibration 
of the beam under end thrusts alone. In this case the expression for y 
becomes excessively large. 

The various terms in equation (10) are now specified except for the con- 
stants in the infinite series. These will be determined by expanding the 
given values of y and dyjdt at ^^ = in a Fourier's series of period tt/Z. It 
appears from the above analysis that the vibrations which occur are three 
in type, viz., those corresponding with the natural period of vibration of the 
beam under end thrust with all the overtones, that corresponding with the 
forced lateral vibration, and that transmitted to it by the periodic bending 
moment at the end. It is to be anticipated on account of this last feature 
that in the case of a spar consisting of any number of bays, the vibrations 
in one bay will be transmitted to any other in the same manner as this 
periodic bending moment operates, and that consequently all the bays vibrate 
with the same periods. 

Effect of Vibrations of Supports.— li in the previous analysis the end con- 
ditions that the two supports are supposed rigid, are modified to correspond 
with the case where these two supports execute vertical harmonic motions 

of the form 

So cos ut and Si cos ut 

respectively, then the analysis will proceed as before, except that condition 
{a) becomes y = So cos iit and Si cos %tt at x = I and a? = respectively. 
From (2) by equating coefficients it will be seen that these imply that the 
term cos ut must arise in (2), and that, therefore, A^ = So and Si at a? = Z and 
a; = respectively, the other values of A^ and B^^ being zero at these points. 
The analysis then proceeds along the same lines and indicates that P„, Qm, E« 
and S„ will be determined as definite finite constants unless the period of 
the vibration of the supports coincides with a free period of the beam under 
end thrust alone, in which case resonance takes place. 

It appears, then, quite generally that resonance, and therefore excessive 
amplitude of vibration of the bay will occur when the period of any one of 
the externally applied forced oscillations become 

2 /ll\ // lol^ \ 

n ' V [gj V \E10iV-^2)r 

where <^^ = FP/EI, and n is any integer. 
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Vibrations of a Wing Sjoar tmder JEJnd Thrust, varying from Bang to Bay 
due to Lift Wires and under Forced Loading due to Wind Forces and Engine 
ThroHbing, — Eegarding the spar as a separate entity in the aeroplane there 
are three principal types of force operating upon it which may set tip and 
affect vibration. The longitudinal thrust, as already explained, is brought 
into being by the lift wires. In addition there are in the longitudinal force 
plane certain bracing wires connected to the spar at intermediate points 
which cause a variation to some extent of the force along each bay. For the 
moment, however, the effect of these longitudinal bracing wires will be 
ignored, and it will be assumed merely that the force along each bay is 
constant, but varies from bay to bay. In normal flight the throbbing of the 
engine sets up a periodic vibration in the structure, and this will be trans- 
mitted by the wires, struts, etc., to the spars. The variation in end thrust 
in each bay due to the effect of this on the wires will be neglected for the 
present, and it will be merely supposed that the result of the throbbing 
will be to impose a periodic displacement upon one or more supports of the 
spar. It is diflScult to estimate up the relative importance of the effect of 
variation in wind forces such as are caused by gusts or oscillations of the 
machine, but it will be assumed that this is equivalent to a periodic lateral 
loading on each bay, and it will be quite sufficient to assume that the period 
of these forces is the same from bay to bay. It is proposed to investigate 
the effect of these forced oscillations upon the bending moments and stresses 
in the spars and upon the conditions under which instability and resonance 
arise. 

Let {r) and (r— 1) be two neighbouring supports of a spar, and as before 
let P of co-ordinates {x, y) be any point in that bay at time t, then by the 
same method as already adopted the equation (of motion) of the rth bay is 

where X^2 _ F^/E^I^ (15) 

and /i^2 -_. lOrjg'Furl.u (16) 

Fy = the longitudinal thrust in the rth bay, 
Eyl^ = the flexural rigidity, 
Wr = weight per unit length assumed constant, 
Wr = lateral force per unit length = Xcos j^-f-Xo. 
The equation (14) has constant coefficients, and therefore the solutions 
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are additive. The boundary conditions that must be satisfied are as 
follows :~ 

(a) At a? = 0, y :=. drf. -=1 D^ cos qt, 

X == Ir, y = dr-i = 1)^-1 cos qt, 

where d is the deflection of a support measured from some datum line fixed 

in space. 

(&) At ^ = M = M,,, 

These bending moments will be obtained by the appropriate moment equa- 
tions to be given later. 

{c) At t = 0, y and dy/dt are given functions of x. 

Generally, the forced vibrations may have frequencies other than those 
corresponding with cos qt, and, moreover, vibrations may occur which are not 
common to both the lateral loading and the vibrations of the supports. The 
above analysis, however, will include these cases if each frequency is taken 
separately, and the resulting moments, etc., added algebraically, since solu- 
tions of the differential equations are additive. When a vibration is not 
common, it could be treated as though present in both the lateral vibration 
and in the vibrations of the supports, except that one of the amplitudes is 
zero. Under these conditions, it will be perfectly general to consider only 
those problems in which every bay and every support receives forced vibra- 
tions of the same period but of different amplitudes. 

Assume y of the form 

y = tAn cos qj + SBj, sin qj, (17) 

where A„ and B„ are functions of x alone, and the two summations may, for 
the present, be either finite or infinite series. 
Inserting this in the equation of motion, we find 

4- S ( ^-~ 4- X/ —f~^ — fjLr%i^Bn ) sin qj = X cob qt + Xq. 
\ dx^ dx^ I 

Hence equating coefficients of cos qj and sin qnt on both sides of this equation, 
it follows that 

S=+V^"-/xrVB„ = (18) 

for all values of n, and 
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except when j„ = q, in which case 

and when §'„ = 0, when 

dx^ doer 

Condition {a) for y will be satisfied if An and B„ are zero at 5? = and 
X =: Ir for all values of qn except f » = ^ when A = D^ and Dr-i at a? = 
and X = Ir respectively. For condition {h) we note that M^, to be determined 
from the appropriate moment equation, must be expansible in the form 

Mr = ^Tn cos qt -{- l^Tn siu qj, (21) 

and differentiating (17) twice with respect to x,(b) reduces to 

UrA/' = Tny etc., (22) 

and E,I3;." = ^V, etc. (23) 

Condition (c) merely states that %An and %Bnqn are given functions of x. 
The solutions of equations (18) and (19) have already been given. They 
are 

An = Cln Sinh knX + I3n cosh JCnX + Ci^ siu knX + ^n COS kn(^, (24) 

where cc^, ^n, ^n\ 0n are obtained from the conditions that 

An = at a:? = and x = Ir 

and — -~ = ~-~ and ^ J^^ at a; = and x = k respectively, 

while the solution of (20) is 

A = asinhA:'(g4-ffcoshya;4-g^'sinfa'fi3'cosfoc+.-p^. 2tSF^ — F\' (^^) 

where ±kn and ±knh and +&' and ±K 

are solutions of ¥-^\r^k^-fir^qn^ = 

and ¥ -\-\r^k^ — fjbr^q^ = respectively. 

Confining attention for the moment to the case of q^ = q, we note that the 
boundary conditions demand 

At i?) = 0, A = T)r and -— = ——- , 
oj = /y, A = Dr-i and -—^ = • .- ■ . 

(Mr JUffi^ 
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Inserting these conditions in (25), we obtain 

1 



a 



/8 



a 



(Jc^ + ¥^)Bmhk'lr 



¥Dr~i + %~^ k^f {ir)-f'' (lr)-/3 cosh ¥ I, 



1 r (r) 



k^ -f ki^ 



jbjrjL 



+ k^Dr-k^^{o)-^"{o) 



r'-r 



(k^ + h"^) sin hi. 



h'^I>r-i 



0-1) 

ET 



k'^^^Qr)-^-^" {Ir)-^' msM, 



r^T 



^' = ^^["£1/^'"^^"*''^^"^"*"^" ^"^ 



where 
At X =^ ^, 



^{x) = 



jL. 



dKjclx = k'a-\-ka ■^"\^' (0) 
_ (r— 1) 1 [ k'lr kU 



/yErly * k^-\-k''^ Lsinh^'4 sin7<;^r. 
(r) 1 i ^^r _ ^^'^r 



+ 



4ErIr ' k^-\-k'^ L tan My tanhF^ 
+ K,(F,g),say, (26) 

where K does not involve the bending moment, (r) and (?^— 1). This may be 
written in the more concise form 



(r - 1) 21,. (F, q) + (r) $8r (F, q) - K, (F,?) = - (J^)^, 



(27) 



where 



3fr (F, q) 



Cf 1. j fy tf> fClff 

lEirlr * lr^{k^'{-k'^) LsinhF^ mikl^ 

0r Or 



I 



1 



aSr (F, ?) = 



E,I, .■ {e,? + Or"') I sinh 8,.' sin ^^ 

L X. \ Or Or 



E4r ' (<9/ + O;^) Ltan ^, tanh O; 
The slope at the same point regarded as on the previous bay is 

+ (»•) JB.+1 (F, ?) + (r+ 1) §tr+i (F, 2)-KV+i (F, ?) = (— j^^^^. (28) 

By equating these two expressions the extended form of the equations of 
three moments for a vibrating spar under an end thrust which varies from 
bay to bay is, writing m^ for (r)^ 

This extended form of the equation of three moments determines the 
coefficient of cosg?f in the expansion (21) for the complete bending moment 
at any support. The remaining coefficients in the expansion for M are 
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determined in a similar manner from equations (18) and (19). These co- 
efficients, however, are brought into existence by the superposed natural 
vibration of the spar under end thrust, and if it be supposed, as is in practice 
quite legitimate, that slight damping is present, then they die out, and there 
are merely left coefficients due to the external oscillations. It is possible, of 
course, that the shape of the wing section might be such that the air forces 
brought into play as the wing vibrates would be of the nature to maintain 
the motion, and thus counterbalance and even over-balance any damping 
effect, but for the moment this state of affairs is not contemplated. The 
coefficients, 51^, 93r, etc., in the extended form of the equation of three 
moments are easily calculated for any bay from the roots of the biquadratic 
in h, while K^ and KV+i are known functions of the end thrust, lateral 
loading and g[. If the end supports of a system occur at the extreme end of 
the spar so that there is no overhang, the bending moments at these points 
are zero, and by commencing at one end and considering successive pairs of 
bays the bending moments at all other supports are determined in exactly 
the same manner as by the ordinary method of three moments. If, however, 
the end supports are not actually at the end of the bays, but overhanging 
pieces occur, from a knowledge of the loading, etc., on these pieces a slight 
modification would require to be introduced. 

In general, these equations will determine finite and definite bending 
moments at each support for a given period 2irlq of the external impressed 
force and of the motion of the supports, and for given values of F along each 
bay. When resonance occurs so that the amplitude of vibration becomes 
excessively large, the bending moments also will become excessive, and the 
resonance condition which determines the resonating period for the whole 
bay then takes the form 

a3i + S52, 9I2 ,0,0 

2I2 , a52 + S3, 5I3 , 

2t3 , 





0, 



0, 



91.-1 

^^n-l + 53n 



= 0. 



(30) 



A reduction formula applicable perfectly generally to this determinant is 
clearly identical in form with that obtained in the previous communication, 
and need not be discussed here. As in the case of the stationary spar treated 
in that paper, a number of theorems arise corresponding with those dealin 
with clamped supports, etc* 
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Condition (30) gives the resonating period for any value of F, while for 
q=^ Q this condition is that required for crippling of the structure non- 
vibrating. This term corresponding to (20a) is always present, and therefore 
for all values of q^ the crippling load is the same as for the non-vibrating 
spar. 

In cases where the effect of the end thrust on the spar is not great, it will 
be true to a first approximation to put \t^ = Fr/Erlr = 0, in which event 
the terms in the previous discussion become considerably simplified. 

The expressions for 21^ and S5r now take the form 

SiirK </>r \Sm (pr Sintl (prl 

^r = Tf'r ■ i- (rArr - Az ) - (32) 

hjr^y, ^r \tanh (pr tan Cpr/ 

where ^r = fr\/(/^rq) = klq^ . ~g|) • (33) 

This formula is at once analogous to the results obtained in the earlier 
paper, for the case of a n on- vibrating spar under end thrust where %r and 
33r are the sums of the corresponding tension and compression terms used in 
that communication, 



^r^r 



S3, = ^-V(/3, + /3/), 



jrJ~r 



(f)r now being of course lr\/{fJ^r9)' 

As a simple illustration of these results, consider the case of two bays 
only of equal fiexural rigidities and of lengths h and h. If there be no end 
thrust the determinantal condition (30) reduces to 

a5i + «2-0, 



6,6, , 



i/__i L_ V-f-^^ ~ ) = 0, 

(^i\tanh0i tan 01/ 0^\tanh02 tan 02' 



1 , 1 _ 1 , 1 

-i~ -^ — — — . — -)- 



tanh 01 tanh 02 tan 0i tan 02 * 
For three bays of lengths 4, h, k the determinantal condition is 

S3i + a32 2I2 

%2 a52+bi 

or (251 + 252)^-212^. 



0, 
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In the case where ii->-0, i.e., the ends of the bay are clamped, 

and the condition becomes 

1 / 1 ___! y^ 1 / 1 . 1 y 

^2^\tanh ^2 tan ^2/ ^2^\sin ^2 sinh c^2/ ' 
i.e., cos ^2 cosh (^2=1, 

i.e., (f>2 = 37r/2 approximately, 

9 Tif //^El 



^.e., 2' 



4 P V \ to 

which determines the resonating period. 

For the case of a bar of length h, fixed at one end and supported at the 
other, 

where I2 -> and therefore S32 -> 0. 
Hence S5i = 0. 

1 1 _ 

t/.e.f — — - - "~" — — Vj 



giving as a first solution 



tanh <^i tan pi 
01 = 3-93 approx., 



i.e., q 1*57^ a/j 



2/J 



The foregoing analysis has assumed that in each bay the end thrust, weight 
per unit length and flexural rigidity are constant. If this assumption be 
removed the general equation of motion can be obtained by an application 
of the same methods and expressed in the form 

82 



^l54) + ^(^^) + -^ = XC0S^4 



9a?2\ Zx?} 3a? \ Zx} g 

where EI, F, and-z^ are arbitrary functions of x. 

A convenient general form may be found for the solution of this equation 
in a form particularly suitable for graphical integration, and at the same time 
the conditions for critical loading, etc., can be easily deduced. The method, 
although quite general, will be applied in the particular case of the crippling 
of a strut of variable flexural rigidity, but under constant end thrust. 

The above equation then takes the form 

dx^'^m ' (34) 



if there is no lateral load, since M = EI 



d?y 



dx? 

VOL. XCV. — A. 2 
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It will be seen that the assumption that no lateral load exists is in reality 
no restriction as regards the present investigation. The condition for 
instability will be a relation that must be satisfied by F in order that the 
bending moments along the strut become infiDite.. . This is in accordance 
with the point of view developed in the earlier communication. 

If the strut be of length /, and pin jointed or simply supported at the 
ends 0? = Oand x =: l^ then it can easily be verified that3 the solution of this 
equation subject to the external conditions due, say, to eccentric loading 

M ™ Mo at ^ = 0, • M = Ml at x - I, 

is M == So-Si + Ss-Sa-f etc., (35) 

where So = A-f AV^;, (36) 

dx\ :=^^ndx, (37) 

where A and A' are constants to be determined from the end conditions. 

The series may be written in the following converging form evidently 
suitable for graphical integration : 



M = A 



+ A' 



*x 



l-F I dx 





0-^^ Jo Jo^^Jo Jo^-^ 



x 



F dx YTT + l^^ dx :prr dx ^w— etc., (38) 

Jo Jo J^-^-^ ' Jo Jo^'^Uo Jo ^^'^ 



where each term may be obtained from the preceding one by two simple 
graphical integrations. 

This may be written 

M== AU4-A'V. 

Inserting the end conditions 

Mo=-AUo4-A'Vo, 

Ml =AUz + A%, 

but evidently Uo = 1 and Vo = from the expressions for IJ and V. Hence 

A = Mo, 

A^=: A-(Mi-MoUO, 

M =: MoU + (Mi-MoUz)X.. 

The bending moment at any point is thus obtained in terms of two con- 
verging series whose terms alternate in sign. Upper and lower limits may, 
therefore, rapidly be set to the sums of these series. 
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For instability, M is to become infinite, and this can only take place when 

Yi = 0, (39) 

unless (Mi-MoUO = 0. (40) 

The conditions for crippling, viz., V^ = 0, is thus 



l-¥ 



Jq -E^J- Jo Jo-*^-*- ^Jo io-*^ 



JjIAiA/ j 



_ V —etc. ... = 0, (41) 

which is perfectly general for any variation in EI, and may be rapidly solved 
for F by graphical integration. 

As a verification, suppose EI constant, then the series becomes 

F P F2 ^^_F^j7 _ 

/ L I - ' ' ' "' *' 1 '",.1*."! ' - •— -** . .""."."""'7'* 'L' 1^ v-/ 'iV^B '^ ' \J * 

EI3! EI5! EI7! 
or multiplying throughout by ^(F/EI) and writing 

<^-S+S-ete- = 0, (42) 

ie., sin (^ = 0, 

Since F =0 is not a solution, 

^y^(F/EI) = ^ = TT, 27r, etc., 

the usual Eulerian result. 

It is clear that the criterion used in the foregoing analysis for crippling, 
namely, that the bending moments become infinite, can be applied directly to 
any case where the solution of equation (34) is completely known. A 
typical case of this nature, well worth analysis in detail, is that afforded 
by the assumption that the variation in the square root of the flexural 
rigidity is parabolic, e,g.y 

EI = EIo f 1 --^T= '^ {e?^f^?f (43) 

At x = 112, I, = I„(l-F/4c^). (44) 

The equation for solution for constant end thrust is 

d?M. , F . c* . M r. /^P-x 



cIq(? EIo (c^— 03^)^ 
This equation is completely integrable, giving 



M = (c2- ^2)1/2 



^/d^Y + B^''"''^^^ 



(46) 



fi — XJ \C-{-Xj 

where X =. JC^-l). (47) 



VIeT 



»-Q' 



2 2 
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Inserting the end conditions, 

M :=z Miat //2') 

M =r Ms at 112 
the condition that M becomes infinite as before takes the form 



; J 



,2c-/, 



■2_c+j\^ 

'2c-lY 



.2c + l/ 



0, 



t,6n 



or 



2c + Ao- //Fc^ i\ 1 



^ 



2t 



2c-l^ 

\Elo 




l)log(g±i)=l. 



(48j 



Hence, solving for F, 



Vi 



>EIo 



1 ) . log t — UTT ; 



2c^/ 



EIo 



[•4 



?^7r 



2- 



log 



2c^l. 






(49) 



Or putting 



E = 4EIote 



nir 



2n 



1 + 



1-f y/q-^r) 



(50) 



l^liQ limit of this expression when r = 1 corresponding with uniform section 
can easily be shown to give 

which is Euler's value. 

The case when 72, ™ must be excluded, for this would give 

and the solution obtained is illegitimate, the real bending moments not 
becoming infinite. The first crippling load is thus 

(51) 



Fr=iM(l-r) fl-f 



TT' 



/ 



2 



1 ,'l + W'(l-r)Y>' 



where r^ =:: Ii /Iq. 

When r = 1/2, so that the moment of inertia at the ends is only 1/4 that 
at the mid-section, this gives FP/¥Ao = 8*28. 
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For a uniform strut whose moment of inertia is only 1/4 that at the mid- 
section of the above strut 



IT 



986. 



EIo 

Hence by tapering the strut to a moment of inertia at the tip equal to 1/4 
that at the central section, the strength has been reduced by 16 per cent. 
The variation in strength does not appear to be greatly affected by moderate 
tapering. The last formula, therefore, should cover most practical cases. 
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Two-dimensional Solutions of Poissons and Laplace's 

Equations, 

Part I. — By Leonard Bairstow, F.E.S., C.B.E. Part II. — By Arthur 
Berry, M.A., Fellow of King's College, Cambridge. 

(Received January 13, 1919.) 

In the course of a preliminary investigation of fluid motion when viscosity 
is taken into account, it was found that a solution of Laplace's equation was 
obtained in very general form. The solution gives exact analytical expres- 
sions where the conditions of the problem admit of exact solution, but the 
chief claim to notice arises from the extremely general character of the 
method if graphical and mechanical integration are employed. 

The problem has been approached from the engineering standpoint, and 
the absence of any proof of the convergency of the series developed will 
immediately be evident to the mathematician. The justification by trial in 
particular cases is all that is attempted at the present moment, but in those 
cases the result is wholly satisfactory. For the analytical solution of the 
problems I am indebted to Mr. Arthur Berry who has added a second part 
to this paper, and has generally acted as helpful critic. 

For single closed boundaries of any shape, or for a combination of separate 
boundaries, the problem of finding a function which satisfies either Poisson's 
or Laplace's equation, admits of direct attack for the region inside a single 
boundary or between boundaries if there be more than one. For the region 
outside a single boundary there is a further condition to be satisfied infinitely 
far from the boundary. 

The problems arising in calculating the torsion of solid or hollow cylinders, 
the flexure of beams, and the motion of an inviscid fluid round a cylinder 



